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1. Introduction 

Recent progress in parallel computing, as well as theoretical advances in the formulation of 
lattice gauge theories with fermions, have allowed the widespread adoption of simulations using 
dynamical light quarks, leading to a significant reduction in systematic errors by removing the 
uncontrollable error inherent in the quenched approximation. 

The Fermilab Lattice, MILC and HPQCD collaborations have an ambitious program which to 
date has made several high-precision predictions from unquenched lattice QCD simulations [1,2]. 
In this, we rely on the Symanzik-improved staggered-quark formalism, specifically the use of the 
asqtad [3] action. While this approach requires the use of the fourth root of the staggered quark 
determinant, all the available evidence is consistent with the conclusion that the resulting theory 
is in the same universality class as continuum QCD, as long as the chiral limit is taken after the 
continuum limit [4]. 

Recent studies of the heavy-quark potential in full QCD [5] have shown an unexpected ap- 
parent increase in scaling violations compared to the quenched approximation. A possible reason 
for this would be that these scaling violations arise from the mismatch between the inclusion of 
sea quark effects in the simulation and the omission of sea quark effects in the improvement co- 
efficients in the action. This mismatch would appear to spoil the (ff{a^) improvement at the level 
of ff{asNfa^). While a systematic study of ff{asa^) effects is generally beyond the scope of the 
current perturbative improvement programme, it is still important to bring up-to-date the calcu- 
lations of the Liischer-Weisz improved gluonic action [6, 7] to include the effects of dynamical 
quarks. This is important also because the Liischer-Weisz improvement is currently included in 
many unquenched simulations [3]. Since the lattice spacing scale is set by measurement of the 
heavy-quark potential, there will be an induced &{asNfa^) artifact by omitting the corrections due 
to unquenching. While such errors are generally smaller than other systematic errors in current 
state-of-the art studies, it is simple to remove them using the result of the perturbative matching 
calculations outlined here. For details, the reader is referred to our paper [8]. 

2. On-shell improvement 

The Liischer-Weisz action is given by [9] 




(2.1) 



where P, R and T are the plaquette, rectangle and "twisted" parallelogram loops, respectively. The 
requirement of obtaining the Yang-Mills action in the continuum limit imposes the constraint 

co + 8ci + 8c2 = l, (2.2) 

which fixes cq given the other two coefficients. This leaves us with ci and C2 to be determined in 
order to eliminate the ff{a^) lattice artifacts. 

If we have two independent quantities Qi and Q2 which, at each order in perturbation theory, 
can be expanded in powers of (/ia), where /i is some energy scale, as 

Qi^ Qi + Wi{iiaf ^dijCj{p.af ^ e {{p.af) , (2.3) 
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then the (ff(a^) matching condition reads 

dijCj = -Wi . (2.4) 

Since this equation is linear, we can decompose the Wi into a gluonic and a fermionic part as 
Wi = wf^^^ ^ NfW^^^^^ and obtain the same decomposition for the c/; thus, especially we do not 
need to repeat the quenched calculation in order to obtain the ff{Nf) contributions.^ 

At tree-level, there are no fermion loops to consider, and hence the tree-level coefficients 
remain unchanged compared to the quenched case [6]: 

ci = - ^ , C2 = 0. (2.5) 



3. Lattice perturbation theory on a twisted lattice 

In lattice perturbation theory, the link variables Uj^ G SU {N) are expressed in terms of the 
gauge field G su{N) as 

C/^(x) = exp (^gaAf^ (^x+ ^fi^^ (3.1) 

which, when expanded in powers of g, leads to a perturbative expansion of the lattice action, from 
which the perturbative vertex functions can be read off. 

As in any perturbative formulation of a gauge theory, gauge fixing and ghost terms appear 
in the Fadeev-Popov Lagrangian; an additional term arises from the Haar measure on the gauge 
group. Here we will not have to concern ourselves with these, since for our purpose we only need 
to consider quark loops. 

To handle the complicated form of the vertices and propagators in lattice perturbation theory, 
we employ a number of automation methods [11,12] that are based on the seminal work of Liischer 
and Weisz [6]. Three independent implementations by different authors have been used in this work 
to ensure against programming errors. 

We work on a four-dimensional Euclidean lattice of length La in the x and 3; directions and 
lengths L^a, L^a in the z and t directions, respectively, where a is the lattice spacing and L.L^^Lt 
are even integers. In the following, we will employ twisted boundary conditions in much the same 
way as in [6, 7]. The twisted boundary conditions we use for gluons and quarks are applied to the 
{x^y) directions and are given by (v = x^y) 

U^{x + Lv) = ayU^{x)a-' , (3.2) 
^(x + Lv) = ay^^{x)^l-^ , (3.3) 

where the quark field ^sc{^) becomes a matrix in smell-colour space [13] by the introduction of 
a "smell" group SU(Ns) with Ns = N in addition to the colour group SU(N). We apply periodic 
boundary conditions in the (z,/) directions. 

^Although doing so provides a useful check on our methods, and we have in fact successfully reproduced the results 
of [7]. 
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These boundary conditions lead to a change in the Fourier expansion of the fields: in the 
twisted {x^y) directions the momentum sums are now over 

NL ^ ^ 

py^mny, <ny<^, v ^ [x,y) , (3.4) 

where ^ = ^ • The modes with (% — fiy — Q mod N) are omitted from the sum in the case of the 
gluons. The momentum sums for quark loops need to be divided by N to remove the redundant 
smell factor. 

The twisted theory can be viewed as a two-dimensional field theory in the (z,r) plane by 
considering the modes in the twisted directions as Kaluza-Klein modes. Denoting n = (nx^riy), 
the stable particles in the (z,/) continuum limit of this effective theory are called the A mesons 
(n = (1,0) or n = (0, 1)) with mass m and the B mesons (n = (1, 1)) with mass V2m [7]. 

4. Small-mass expansions 

To extract the (ff{a^) lattice artifacts, we first expand some observable quantity Q in powers of 
ma at fixed ruqa: 

Q{ma,mqa) = a^^\mqa) +a^^\mqa){ma)^ + ^ ((ma)^, (ma) log (ma)) (4.1) 

where the coefficients in the expansion are all functions of rUqa. There is no term at (ff [{maY log (ma)) 
since the gluon action is improved at tree-level to 0{a^) [7]. Although we ultimately wish to ex- 
trapolate to the chiral limit, we cannot set m^a = straight away, since the correct chiral limit is 
rUqa ^0, ma ^ 0, m^/m > C, where m = ^ as before and C is a constant determined by the 
requirement that a Wick rotation can be performed without encountering a pinch singularity. This 
requires us to consider a double expansion in mqa.ma and carry out the extrapolation to m^a = 
for the coefficients in Eqn. (4.1). 

To extrapolate to the chiral limit, m^a 0, we will fit the coefficients in the expansion for Q 
in ma to their most general expansion in m^a for small m^a. 

For ^Q^^ (m^a) we have 

af^ {mqa) = hf^ log(m,a) + af^ . (4.2) 

Since we expect a well-defined continuum limit, d>^^ (mqa) cannot contain any negative powers of 
mqa but, depending on the quantity Q, it may contain logarithms; b^^^ is the anomalous dimension 
associated with Q, and can be determined by a continuum calculation. 
For ^2^^ (m^a) we find 

(mqa) = + 4^ + US + 45^ log(m^a)) (mqaf + ^ ((m^a)^) . (4.3) 

After multiplication by (ma)^ the {mqa)~^ contribution gives rise to a continuum contribution to Q, 
and a^^}_2 calculable in continuum perturbation theory. There can be no term in {mqa)~^ log(m^a) 
since this would be a volume-dependent further contribution to the anomalous dimension of Q, and 
there can be no term in \og{mqa) since the action is tree-level 0{a^) improved. 

In the chiral limit m^ 0, the term wi that appears on the right-hand side of Eqn. (2.4) is a^^^, 
and it is this limit and this coefficient that we will concern ourselves with hereafter. 
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(a) 



(b) 



Figure 1: (a) The fermionic one-loop diagrams contributing to the A meson mass renormalisation as well 
as to the wavefunction renormalisation for A and B mesons, (b) The fermionic diagrams contributing to the 
irreducible three-point function. 

5. The A meson mass 



The simplest spectral quantity that can be chosen within the framework of the twisted boundary 
conditions outlined above is the (renormalised) mass of the A meson. The one-loop correction the 
the A meson mass (for A mesons with positive spin) is given by 



m 



(i)_ 



-Zo(k) 



2m 



(0) 



(5.1) 



k^iim^^\0,m,0) 



where Zo(k) = 1 + ((ma)^) is the residue of the pole of the tree-level gluon propagator at spatial 
momentum k, and mf"^ is defined so that the momentum k is on-shell. We consider the dimension- 
less quantity m"^^ /m. The fermionic diagrams that contribute to this quantity are shown in figure 1 
(a). 



The anomalous dimension of ma is zero and so using Eqn. (4.2) we have b\ 



K,i) 

0,0 



0. From 



gauge invariance we find a2^^2^ = and a^^^'^^ {^qo) = 0, which together with the previous result 
implies that a^^^'^^ = 0. 

The & {as{mdf') contribution from improvement of the action is given by [7] 



m 



(1) 



Aimp^ = -(4^^ -4^^)(ma)^ + ^((ma)4) . 



(5.2) 



6. The three-point coupling 



An effective coupling constant A for an AAB meson vertex is defined as 

A = ^ov/Z(k)Z(p)Z(q)e,ri'2'^(/:,p,^) 



(6.1) 
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where we have factored out a twist factor of ^Tr([ry^,rp]r^) from both sides, and the momenta 
and polarisations of the incoming particles are 

k={iE{k),k) p={-iE{p),p) q={0,q) e= (0,1,-1,0) ^^^^ 

k = (0,m,/r) p = (m,0,/r) q = (— m, — m, — 2/r) 

Here r > is defined such that £'(q) = 0. This coupling is a spectral quantity since it can be related 
to the scattering amplitude of A mesons [10]. We expand Eqn. (6.1) perturbatively to one-loop 
order and find (up to ff{{maY) corrections) 



Ad) 1 .\r(i) 4 d 



1 A d 



-[^-i-X]- ^ir<i\k) - ( 1 - i-m^ ^ ( e'e^n^\q) 

m \ 24 I m kndkn ' \ n ) ^^2 \ ij vy; 



72 



(6.3) 



kQ^iE{^) \ \2 J dql V 'J 

The fermionic diagrams contributing to the irreducible three-point function T^^^ are shown in figure 
1 (b). Continuum calculations of the anomalous dimension and infrared divergence give 



b^^^'^ - -^2^ a^^^'^ - -^^2" (6 4) 



The improvement contribution to X is [7] 



A- — =4(9cS^^-74^^)M^ + ^((ma)4) . (6.5) 



"imp 



m 



7. Continuing to imaginary momenta 

The external lines of the diagrams are on-shell, but with complex three-momentum k; in the 
Euclidean formulation is also imaginary. In evaluating the loop integrals that are not pure tad- 
poles, care must be taken to ensure that the amplitudes calculated are the correct analytic continu- 
ations from the Minkowski space on-shell amplitudes defined with real three-momenta to the ones 
in Eqn. (6.2). 

The situation is complicated by the presence of two mass scales m,m^. The integrals are 
evaluated after performing a Wick rotation in k{), taking care to avoid contour crossing of any poles 
that move as r is continued from r = to r = m/\/2. This requires mq/m>C, where C is a constant 
dependent on the graph being considered. After the Wick rotation in the (Euclidean) integration 
contour for (o^:, in one case, ^3) must be shifted by an imaginary constant. 

8. Results 

To extract the improvement coefficients from our diagrammatic calculations, we compute the 
diagrams for a number of different values of both L and niq with Nf = I, N = 3. At each value 
of ruq, we then perform a fit in ma of the form given in Eqn. (4.1) to extract the coefficients 

a^'^\mqa), ^ = 0,2. 

Performing a fit of the form (4.2) and (4.3), respectively, on these coefficients, we get the 
required coefficients of the (ff{a^) lattice artifacts in the chiral limit to be 

4"^^'^) = 0.00361(1) (8.1) 



^2,0 
^2,0 



-0.140(1) (8.2) 
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Figure 2: (a) A plot of the fermionic contributions to the one-loop A meson self-energy m>^ jm against 
{md)^ . The vanishing of m>^ jm in the infinite- volume limit can be seen clearly, (b) A plot of a>^'^^ against 
niqa which shows the agreement between the numerical lattice results and the known anomalous dimension, 
(c) A plot of ^2 against m^a with the analytical continuum result for the infrared divergence shown for 
comparison. 



These coefficients are to be identified with the of Eqn. (2.4). 

Solving equation (2.4) for c^^\ our results can be summarised as 

c^l^ = -0.025218(4) +0.00486(13)7V/ 
-0.004418(4) +0.00126(13)/// 



(8.3) 
(8.4) 



where the quenched (A// = 0) results are taken from [7]. With A^/ = 3 the shift from the quenched 
values is surprisingly large, and may have a significant impact. 
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